I.V. Stasyuk, V.O. Krasnov 



o 

(N 



3 



O 

o 



> 
o 

O 

m 



x 

S3 



PACS 71.10.Fd; 71.38 



I.V. STASYUK, V.O. KRASNOV 

Institute for Condensed Matter Physics of the Nat. Acad, of Sci. of Ukraine 
(1, Svientsitskii Str., Lviv 79011, Ukraine) 

ENERGY SPECTRUM OF THE PSEUDOSPIN-ELECTRON 
MODEL IN A DYNAMICAL MEAN-FIELD APPROACH 



The pseudo spin- electron model in the case of infinite on-site electron repulsion is investi- 
gated. The electron energy spectrum is calculated within the framework of the dynamical mean 
field theory (DMFT), and the alloy analogy approximation is developed. The effect of the 
pseudo spin- electron interaction, local asymmetry field, and tunneling-hke level splitting on the 
existence and the number of electron subbands is investigated. The relation of the pseudospin- 
electron model to the problem of energy spectrum of boson-fermion mixtures in optical lattices 
is discussed. 
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1. Introduction 

The pseudospin-electron model (PEM) is one of the 
models, which are used in the physics of strongly cor- 
related electron systems in recent years. Application 
of the model to high-temperature superconductors al- 
lows one, for example, to describe the thermodynam- 
ics of an anharmonic oxygen ion subsystem and to 
explain the occurrence of inhomogeneous states and 
the bistability phenomena see also @). In this 
model, one considers the dynamics of locally anhar- 
monic structure elements (using the pseudospin vari- 
ables to describe them), interaction between pseu- 
dospins and electrons, and the asymmetry of local 
anharmonic potential wells. The electron subsystem 
is described by the Hubbard Hamiltonian. 

In Q , a pseudospin-electron model of ion intercala- 
tion in crystals was formulated; the pseudospin repre- 
sentation is also used in the description of ionic con- 
ductors • The thermodynamics of the model has 
been investigated in the mean-field approximation; it 
was shown that a new phase with (S x ) ^ appears 
due to the ion hopping between sites. This new phase 
is an analogue of the superfluid phase in the systems 
of hard-core bosons or the superionic phase in ionic 
(protonic) conductors. 

The pseudospin-electron model in its general form 
has also a direct relation to the ultra-cold boson- 
fermion mixtures in optical lattices, which are the 
object of an intense theoretical investigation during 
last years In the hard-core boson limit (that 

corresponds to the large on-site repulsion between 
Bose atoms), when the pseudospin representation is 



© I.V. STASYUK, V.O. KRASNOV, 2013 

68 



used, one has to do with a pseudospin-fermion ana- 
logue of PEM . The transverse field exists here in 
the phase with a Bose-Einstein (BE) condensate and 
is proportional to the order parameter (S x ), while 
the longitudinal field plays the role of a chemical po- 
tential of bosons. In this case, the investigation of 
the electron spectrum within PEM gives, at the same 
time, an information about the fermion spectrum of 
the boson-fermion system with < I in the men- 
tioned phase. 

Investigations of the thermodynamics and the 
dynamics of PEM were performed mostly within 
the Hartree-Fock-type approximation (in the weak 
pseudospin-electron coupling case) or the general- 
ized random phase approximation (at the strong cou- 
pling); the problem is reviewed and discussed in 
[2]. The conditions of the appearance of modulated 
phases, phase transitions between different uniform 
phases, and phase separation were established. 

In the DMFT approach, the Hamiltonian with 
strong correlations is taken in the infinite space di- 
mension (d — > oo) limit; this leads to a reformula- 
tion of the problem and the transition to the solution 
of the single-site problem described by an effective 
Hamiltonian [12hI4 |. Only for simplest cases such as 
mobile particles in the Falicov-Kimball model, one 
can solve analytically this problem. An exact so- 
lution exists also for the pseudospin-electron model 
without transverse field [15|. There are also some ap- 
proximate analytical approaches such as: Hubbard-I, 
Hubbard-III, alloy analo gy ( AA), modified alloy anal- 
ogy (MAA) etc, see (TaLLZl- 

In this work, the alloy analogy approximation for 
the single-site problem is used within the pseudospin- 
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electron model. Our problem is solved in the limit 
of infinite value of single-site electron interaction. 
The previous consideration of this problem in the 
Hubbard-I approximation 18( revealed a complicated 
structure of the spectrum and the presence of some 
number of subbands. In [19| . the electron energy 
spectrum of the pseudospin-electron model allowing 
the interaction of near-energy subbands was consid- 
ered. The effective single-site problem was solved 
within the auxiliary fermion field approach with the 
help of the procedure of different-time decoupling of 
higher order Green's functions [lf|. 

Our main task is to investigate a reconstruction of 
the electron energy spectrum and to describe the ef- 
fect of band splitting at a change of the longitudinal 
and transverse fields and the pseudospin-electron in- 
teraction constant. 

2. Hamiltonian of the Model and its 
Transformation 

The Hamiltonian of the pseudospin-electron model is 
H = Y J H % + tijcttj*, (!) 

<i,j>,cr 



where f2 in the single-site part of the Hamiltonian is 
the tunneling-like level splitting, g is the pseudospin- 
electron interaction constant, and h is the asymmetry 
of the local anharmonic potential: 

Hi = -/i(n ijt +n a )+5(n iit +n iji )5f-OS , f-/i5f. (2) 

The second term in (1) describes the electron site-to- 
site hopping (between nearest neighbours). 

The term U n^n^ is not included, as we investi- 
gate our problem in the limit of the infinite poten- 
tial of single-site electron interaction U — > oo, where 
all states with double cite occupation are absent. 
The single-site Hamiltonian is considered as the zero- 
order one with respect to the electron transfer. It 
is useful to introduce the standard single-site basis 
\R) = \rii f, Tin, S'), with six eigenvectors [Tl |: 



0,0,1 



13} 



o,i, i 



14} 



i,o, i 



i = 



i 



o,o,-- 



13} 



0,1,-: 



1 



14} = 



1,0,-- 



(3) 



Using Hubbard X-operators that act in the space 
of such eigenvectors, we can write down the electron 
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annihilation (creation) operators and the pseudospin 
operators (in the limit, when the state |2) = |1, 1, ±|} 
is unoccupied) as follows: [l8[ 



— Yl4 i vl4 

— A-i + A i , 



(4) 



c^Xf+X 13 

Then, the single-site part of the Hamiltonian can 
be expressed by means of X-operators in the following 
way: 

^ = (-/i+f ~)(Xf + Xf)+ 



2 2' 

-f ) + J(xF-x 11 )- 



+(-M-f + ^)(Xf + X?, 



+ ^(xf + XY+ Xf + Xf + Xf + xf ). (5) 

This Hamiltonian is diagonal in the case O = 0. But 
if the tunneling splitting is non-zero, we have to use 
a transformation 



fR\ _ ( cos0 r sin0 r \ (r\ 
\RJ \— sin0 r cos0 r / \r J 

to diagonalize it. Here, 

, . n r g — h 

cos(20 r ) = — , 

ni = 0, ns = m = 1. 

In that way, we have 

H = Y j e r xr+ £ t ijC fc ja , 

i,r <z, i 7>,(T 



(6) 



(7) 



-1,1 



i,r <i,j>,a 



£ 3,3 ~ £ 4,4 

Here, making transformation (6), we obtain 

= COS 04i (X, 



(8) 



l 4 + Xt)+sm^ 1 (xf-X}% 



= cos 3i (X, 13 + Xf) + sin 031 (X 13 - X 13 ), 



COS 041 = COS (04 — 0i ), COS 03 

sin 04i = sin (04 - 0i), 
where X-operators act on the new basis 



1 = COS (03 - 0i), 

iin0 3 i = sin (03 - 0i). (9) 
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3. Dynamical Mean Field Theory Approach 

The transition to the d = oo limit in the DMFT ap- 
proach is accompanied by the scaling of the electron 
transfer parameter: 



U 



\fd" 



(10) 



In particular, the self-energy part of electron Green's 
function becomes purely local fl3, 14|: 



S>j,<r(w) = T, a 8 



•J • 



The Fourier-transform Y,ij. a (u)) 
momentum-independent: 



(11) 

is, hence, 



S CT (k,o;) = Y. a (w). (12) 
Electron Green's function in the (fc,w) representation 

Gi^^^-^GaA") (13) 



can be expressed as 



tk ' 



(14) 



where S CT (o;) is the part, which is irreducible (in the 
diagrammatic representation) according to Larkin. 
To calculate the E a (uj) function, the effective single- 
site problem is used. The transition to this problem 
corresponds to the replacement 



Texp^ 



o 

JdT'Y,UT-T)at{T)a a {T^ 



drx 



*W), 





where 
Hq = Hi 



(15) 



(16) 



and J a (t — t ) is an effective time-dependent field (co- 
herent potential) that is determined self-consistently 
from the condition that the same self-energy part 
S CT (u;) determines the lattice function G%(u>), as well 

as the Green's function (oj) of the effective single- 
site problem: 
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[E a (uj)}-i-j a (ujy 



(17) 



In this case, we have 

=GuA") = iY, G U")- 



(18) 



The system of simultaneous equations (14), (17), and 
(18) becomes closed, when it is supplemented with 
the functional dependence 



GP(w) = f([ML>)]), 



(19) 



which is obtained as a result of solving the effec- 
tive single-site problem with the statistical operator 
exp(-0H eS ), 

4. Reformulation of Wick's Theorem 
for the Single-Site Problem 

To find out relation (19), let us calculate electron 
Green's function using an expansion in powers of the 
coherent potential J a {ui). In the zero approximation: 

-{TX™{t)X™(t')) = -g q p {T - t')(X«« + XPP) . 

(20) 

In the frequency representation: gQ P (u; n ) = —{iui n — 

^pq) i ^pq ^p ^-q- 

Using Wick's theorem for the Hubbard operators 
12011, we can see that 



X 41 (T')X i4 (r) = -g^( T - T')(X LL + X 44 ) T 



X 41 (r')X 14 (r) = 



X 41 (r')X i4 (r) = -g^(r - r')X 44 (r') 



X^(r')X^(r) = -g^(T-T')X^(r') 

As a result of such a procedure, the Bose-type X- 
operators appear. 

The alloy analogy approximation (see [Til , 17 1) 
means the neglect of all non-diagonal Hubbard op- 
erators in Wick's pairings. Such an approximation 
leads to the next result: 



c+(r')X 14 (r) - - 3o 14 (r - r')(X 11 + X A %, cos<M, 



c+(t')A^ 4 (t)=0, 



c+(t')X 14 (t) = -gl\T - r')(X n + X^) T , cos0 4 i, 
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C +(r')X Ti (r)=0, 



c+(t')X 14 (t) = 3o 14 (r - r')(X n + X A \, sin^i, 



c+(r')X 14 (r)=0, 



c+(t')X 14 (t) = -gl 4 {T - t')(X u + X A %, sin. 



c+(r')X 14 (r)=0. 
Using (9), we have 



(21) 



C +(r')c T (r) = -g^(r - r'){X^ + X A %, cos 2 <£ 41 



-^ 4 (T-r')(X n + X 44 ) r , cos 2 ^41 



-^ 4 (r-r')(X u +X 44 ) r ,siii 2 ^i- 



-9l\r - r')(A n + X AA ) T , sin 2 4 1 . (22) 



This result shows us that, in the case of the alloy 
analogy approximation the pairing of Fermi-operators 
decomposes into the sum of terms that are the pro- 
jections on single-site states (because of the action of 
X rr operators). This is the main difference from the 
case of ideal fermions, where we have Green's func- 
tions as a result of the pairing. 

Now, we can rewrite 



(r')c t (r) = 

-[.9o 4 (r-0 cos 2 041+5o 2 (t-t') sin^X 11 ^')- 
-{gf(T-T>) cos 2 041 +gl\T-T')( S mfa 1 ]X ri (T>)- 
-[glHr-r 1 ) cos 2 04i+5<iV-T') sin04i]* 4 V)- 
-bF(r-r') cos 2 04i+5o 24 (r-r') sin^ijA^V) = 
- 9 $X^(t>) g$X r \r<) - $X*(/)- 
-g®X&{T>) = -Y,$X rr {T')- (23) 

r 
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5. Single-Site Green's Function 

In general, electron Green's function of the effective 
one-site problem reads 

T( _ m (Tc a (r)cj(r')e-^) 



(Tc a {T)ci(T>)a{P)) a 
(*G8)>o 



(24) 



The numerator and denominator in this expression 
will be calculated separately using an expansion in 
terms of the coherent potential J a (r— t'). As the first 
step, we illustrate the second order in this expansion 
with four operators of creation and annihilation of 
electrons: 

(Tc t (r)c+(r / )4(r 1 )c t (r 2 )}o = 



(Tc+(t)c t (t)c+(ti)c t (t2)) + 



+(T C |(r 1 )c t (r) C {(r')c t (T 2 ))„ = 

= -5>#(r - t') 9 $(t 2 - r 1 ){X rr )o+ 

r 

+ E»of ( T " ri)g$fo - r')(X rr )o (25) 



and 



(r CT (r)4(r') C +(r 1 ) C ;(r 2 ))o 



= -(T4(Ti)c t (T)c+(T 1 )c l (T 2 )) Q = 
= -E-9ot (t Offfffa - n)(X"V (26) 

r 

Here, the pairing of Fermi-operators is performed ac- 
cording to (25). The diagonal A-operators, which 
appear during this procedure, we multiply, by using 
the rule X rr XPP = X rr S rp . As a result, only the 
averages (X rr )o are present. 

We can also consider the third order in our expan- 
sion with six operators of creation and annihilation 
of electrons: 

( rc t( T ) c t ( T ') c t ( T i)c t (T 2 )4 (T 3 )c t (T 4 ))o = 



I ■ 1 

= -< Tc t ( T ') c t( T ) c t (n)c t (T 2 )c| (r 3 )c t (T4))o+ 



+(rc+(T')c+(rl)c t (r) CT (T 2 )c+(r 3 )c t (r 4 ))o+ 
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+ (T C +(r')c+(T 1 )c t (T 2 )) oC +(T 3 )c t (r)c T (r 4 )) 



-(Tc+(T') C t(r)c+(T 1 ) CT (T 2 ) C +(r 3 )c t (T4))o+ 



i — * — i i ■ — i i — 



+ (TcJ (r')c t {T)cJ (ti)c t (t 4 )c^(t 3 )c t (t 2 ))o+ 



Ct(T2)4(Tl)c t (r)c+(T 3 )c t (T 4 ))o- 



I m 1 



(Tc+(T')c t (T 4 )c+(T 1 )c t (T)c+(T 3 )c t (r 2 )) + 



—\ i 1 r~ 

+ ( Tc t ( r ') c t(T4)4 ( t i)c T (t 2 )4 (r 3 )c t (r))o- 



I * 1 



-(Tc+iT'^iT^C+iT^C^C+iT^T^o. (27) 

Finally, 

(Tc t (r)4(T / )4(Ti)c t (T2)c|(T3)c t (r 4 ))o = 

= £>#(t - t')#(t 2 - njfffffa - r 3 )(X")o- 

r 

- $>S } (t - r')4 ) (r4 " r 1 )4 ) (r 2 - r 3 )(X-) - 

r 

- $>o? (r 2 - r')#(r - n^fa - r 3 )(X"')o + 

r 
r 

E# fa - rOflff^ - n)#(r - T 3 )(X"') + 

r 

+ E-9o , tV 2 -r')5^(T-r 3 ) ff W(r 4 -r 1 )(X-) . (28) 

r 

The similar procedure is also actual in the case of 
higher order terms. Using the diagrammatic series, 
we can separate the connected and disconnected "vac- 
uum" (without external vertices) parts of diagrams. 
The former form a geometric progression in the fre- 
quency representation. The latter look like closed 
rings of different lengths (created by unperturbed 
Green's function and coherent potential lines) and 
give exponential contributions in subspaces \r) after 
the summation of infinite series. 

So, the numerator in Green's function 
(Tcf(T)c+{T')) nunl reads 

(Tci-(T)c^(T')) num 
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+ffot ( Un ) J t ( w ™ ) 9 of {u n )J t {u! n ) g$ (u n ) - 



~? 1 + .9ot( w «) J t( w «) 



(X^)oe^. 
(29) 



Here, Q r in the analytical form is 
Qr = ^y^2,g { ^{uj n )J <7 (uj n )- 

~ 2 [ EE Soa ( w ™ ) Jct ( w ™ ) 
+ 3 [ E E 9oJ ( w « ) J <? ( w ™ ) 

= EE lll (i + .9o r , ) K)^K))- 

In particular, 

= E ln ( x + .# T) K) J t K))+ 

+ ^ln(l + .g^ 1) K)J i ( W „)), 

^3,3 = E ln ( x + k)^K)), 

<3 4 .4 = E ^ + ff0f ( W »») J t( W n))- 

The next step is to calculate the denominator 

(5(/3))o = l - J dm / dT 2 E J -fa - T 2) x 



(30) 



(31) 







x(t c +(t 1 ) Cct (t 2 )) + - y dn... 



X EE Jff ( Tl ~~ T 2) J ff'( T 3 - T 4)X 

x(T C +(t 1 ) Cct (t 2 ) C +(t 3 ) C ^(t 4 ))o- 



dr 4 x 



(32) 



In the diagrammatic representation, this series is ex- 
pressed through the set of "vacuum" diagrams. The 
final result could be expressed in terms of contribu- 
tions Q r of the above-mentioned ring diagrams. 
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1,5 
CO 
1.0- 



o.o 




Im J(w) > 



41 
41 



h 1.5 



Fig. 1. Dependence of electron band boundaries on the asym- 
metry of the local anharmonic potential h (g = 0.5, O = 0.1, 
T = 0.02, fj, = 0, W = 0.5). Hereinafter, dashed lines repre- 
sent the energies of the transitions between single-site electron 
levelspq = \ pq = e p — e q without hopping 

In such a case, we have 



r(/3)) = l + J2[Qr + l i Q 2 r + l i Q 3 r - 
r 

r 

Finally, our analytical result is 



(33) 



(Tc 



where a =f or \.. 



l+ffiy(w»)Jo-(Wn) 



p 



(34) 



6. Electron Energy Spectrum 

Now, we have the closed system of equations to cal- 
culate Green's function G^' (lu) and the coherent po- 
tential J-f(u>): 



G^{w) 



[S CT ( w )]-i - J a (cj) N 



k 



E 9 0a ("nj / yrr\ B 



p 

Q r = EE ln ( X +flte(w«)^K)). (35) 
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Im J(o>)=0 



,41 

*— — 

41 
41 



Im J(m)=0 




0,0 0,5 1,0 h W 

Fig. 2. Dependence of electron band boundaries on the 
asymmetry of the local anharmonic potential h (g = 0.5, Q = 
0.3, T = 0.02, n = 0, W = 0.5) 



CO 

1,0 

0,5- 



lm J(a>)=0 



41 



41 



Im 



J(b>)=0 



41 



»41 



Fig. 3. Dependence of electron band boundaries on the 
asymmetry of the local anharmonic potential h (g = 1.0, Q = 
0.1, T = 0.02, fi = 0, W = 0.5) 



Here, to calculate (X rr )o, we use the iterative pro- 
cess: the next (X rr )o value depends on the previous 
(X rr } one as 



(X rr ) 



(X^) e c 



J2 jXPP) eQ P ' 

p=l,l,4,4 



The initial averages are taken from the Boltzmann 

distribution (X rr ) = ^ e -7e„ ■ 
p 

To sum over k, we use the semielliptic density of 
states po(t) = -^^i\/W 2 — t 2 . In this case, J a {uj) = 

^-G^ito) [131; and our final equation for the coher- 
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-1,2 -| 1 1 1 1 1 1 1 1 

■0,5 0,0 05 1,0 h 1,5 



Fig. 4- Dependence of electron band boundaries on the 
asymmetry of the local anharmonic potential h (g = 1.0, Q = 
0.3, T = 0.02, n = 0, W = 0.5) 




1 1 1 [ 1 ] ' 1 ' 1 ' 1 ■ r~ 

0,0 0,2 0,4 0,8 0,8 1,0 1,2 1,4 



Fig. 5. Dependence of electron band boundaries on the 
pseudospin-electron interaction constant g (h = 0.1, Q = 0.1, 
T = 0.02, n = 0, W = 0.5) 




-! . 1 1 1 1 [ 1 1- 

0,5 1.0 1.5 2,0 ri 2,5 



Fig. 6. Dependence of electron band boundaries on the 
pseudospin-electron interaction constant g (h = 1.0, fi = 0.3, 
T = 0.02, n = 0, W = 0.5) 




Fig. 7. Dependence of electron band boundaries on the 
tunneling-like level splitting U (g = 1.0, h = 0.1, T = 0.02, 
M = 0, W = 0.5) 



ent potential J a {uj) is as follows: 



W 2 ^ l+fl& ) (o>»)-M<"») 



rr\ r „Q r 



(X rr } 



p 



(36) 



In a usual way, we perform the analytical continu- 
ation on the real axis (iu> n — > uj — i5), and only the 
solutions with $tJ a (u)) > must be considered. 

Electron band boundaries are determined from the 
condition QJ a (uj) — > 0. Their dependences on the 
asymmetry of the local anharmonic potential are 
shown in Figs. 1-4. 

One can see the effect of the tunneling-like level 
splitting Q on the width of existing bands (we can 



compare the plots with constant g = 0.5 and to see 
that the increase of f2 from 0.1 to 0.3 leads to a not- 
icable broadening of bands) . 

One can see also how the pseudospin-electron in- 
teraction g leads to the appearance of two additional 
bands, when h ~ f, and one band otherwise (this 
can be seen with comparing the plots with different 
<7, but with the same £1). The dependences of elec- 
tron band boundaries on the pseudospin-electron in- 
teraction constant g are shown in Figs. 5-6. One can 
see that there exists some critical g ~ W, at which 
the additional gap appears in the spectrum. This 
result is similar to the Mott transition type recon- 
struction obtained in (l5l | for the pseudospin-electron 
model without tunneling-like level splitting. 
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lm J(m) 



!m J(a 



g=i.o 






lm J(co) 



g=L5 





-1,0 -0,5 0,0 0,5 



1,0 CO 1,5 



Fig. 8. Electron density of states at different values of the pseudospin-electron interaction constant g (T = 0.02, /i = 0, W = 0.5) 



Similarly to the dependence of the local anhar- 
monic potential on the asymmetry, we can see also 
how an increase of leads to the band width growth. 
At g rj 2h, the additional bands appear in the 
spectrum. 

The dependences of electron band boundaries on 
the tunneling-like level splitting Q are shown in Fig. 7. 
Here, one can see the presence of two critical val- 
ues w 0.3 (band splitting) and il w 0.6 (band 
emerging) . 

The densities of states p a {uj) = ^ImG^ w — i0 + ) 
for different values of the pseudospin-electron inter- 
action constant g are shown in Fig. 8. Initially, we 
observe the splitting of one band at g « W. After- 
wards at g w 2h, two central bands emerge, and the 
further increase of g gives four bands in the electron 
spectrum. 

7. Conclusions 

1) The electron energy spectrum of the pseudospin- 
electron model is considered. For this purpose, the 
dynamical mean field method is applied. The effective 
single-site problem is solved within an original ap- 



proach based on the use of a generalization of Wick's 
theorem. The alloy analogy approximation is used. 

2) Electron transitions with a possibility of the si- 
multaneous flip-over of the pseudospin (or at its un- 
changed orientation) manifest themselves in a com- 
plication of the electron spectrum which consists in 
the appearance of additional subbands. The effects 
of bond splitting and creation of new gaps take place, 
depending on the longitudinal field h, interaction con- 
stant g, and transverse field values (Figs. 1-7). 

3) The rearrangement of the spectrum, similarly to 
the Mott transition, occurs not only due to the on-site 
peudospin-electron interaction, but also under influ- 
ence of the transverse field that intertangles the states 
with different pseudospin orientations. 

4) The obtained results show that, in the real sys- 
tem, which exhibits the local anharmonicity of lattice 
vibrations, the metal-insulator transition due to the 
short-ranged electron correlation is influenced by the 
anharmonic subsystem. Changing the parameters of 
local anharmonicity (e.g., the potential well shape), 
one can affect the conditions of the appearance of a 
a gap. 
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5) The obtained results (Fig. 7) point out, in par- 
ticular, that the additional gaps can appear in the 
fermion band spectrum at the increase of f2. We have 
considered the case where, at Q = 0, the spectrum is 
already split. The problem in such a limit is reduced 
to the Falikov-Kimball model with a Mott gap at the 
chosen values of model parameters. In the case where 
PEM is applied to the boson-fermion mixtures on a 
lattice, the increase of the f2 parameter corresponds 
to the deepening into the phase with BE-condensate 
(il ~ Yl^ij(^j))- Additional gaps could appear at 
certain critical values of the order parameter ( S x ) . 
This problem (which now attracts an attention [2l| ) 
requires, however, a more detailed investigation. In 
the present study, based on the standard version of 
PEM, the boson transfer is taken into account in the 
mean field approximation. At the same time, the con- 
tributions of the collective pseudospin dynamics to 
the fermion spectrum can be essential. 

6) It should be mentioned that we use an ap- 
proach based on the formalism of X-operators. 
Such a scheme gives a possibility to consider, in a 
unified way, the creation of composite excitations 
(fermion+boson; fermion+boson hole) and their con- 
tributions to the total spectral density. 

7) The more complete analysis of a reconstruction of 
the energy spectrum will be a subject of our subse- 
quent consideration. It is referred, in particular, to 
the region of half- filling (// ~ 0, h ~ |), where the 
instability with respect to the appearance of a mod- 
ulated (CDW-like) phase takes place in PEM at low 
enough temperatures. In such a case, the additional 
complication of the spectrum will arise as a result. 
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EHEPrETHHHHH CI1EKTP 

nCEB/lOCniH-EJIEKTPOHHOI MO/1EJII 

B METO^I AHHAMIHHOrO CEPE^HbOrO IIOJI5I 

P e 3 kj m e 

/locjiifl^KeHO nceBflOcnm-ejieKTpoHHy Moment y BHnaflKy He- 
CKiHMeHHo'i B3aeiviOflii BifliHTOBxyBaHHH ejieKTpomB Ha By- 
3jii. EjieKTpoHHHii eHepreTHHHHH cneKTp MO^ejii po3paxoBa- 
ho b paMKax MeTO/ry r a, HHaM i HHOro cepe^HBoro nojin (^CIl) 
Ta Ha6jiH>KeHH5i cnjiaBy. ^UpcjiifljKeHO bitjihb nceB^ocnm- 
ejieKTpoHHOi B3aeiviOflii, jiOKajiBHoro nojin acHMeTpi'i Ta Ty- 
HCjiBHoro po3inenjieHHH piBHiB Ha icHyBaHHH Ta kIjibkIctb 
ejieKTpoHHHx ni^soH. 06roBopeHO mo^kjiubIctb 3acTOcyBaHH5i 
nceBflOcnm-ejieKTpoHHoi - MO/jejii flo po3paxyHKy eHepreTHHHO- 
ro cneKTpa 6o30H-d>epMioHHHx cyMimepi b ohthmhhx rpaTKax. 
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